Quantum Superinductor with Tunable Non-Linearity 



(N 

o 

oo 



M.T. Bell, 1 LA. Sadovskyy, 1 L.B. Ioffc, 1 ' 2 A.Yu. Kitaev, 3 and M.E. Gershenson 1 
1 Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08854, USA 
2 LP THE, CNRS UMR 7589, 4 place Jussieu, 75252 Pans, France and 
3 Caltech, Institute for Quantum Information, Pasadena, California 91125, USA 

We report on the realization of a superinductor, a dissipationless element whose microwave 
impedance greatly exceeds the resistance quantum Rq. The design of the superinductor, imple- 
mented as a ladder of nanoscale Josephson junctions, enables tuning of the inductance and its 
nonlinearity by a weak magnetic field. The Rabi decay time of the superinductor-based qubit ex- 
ceeds 1 ps. The high kinetic inductance and strong nonlinearity offer new types of functionality, 
including the development of qubits protected from both flux and charge noises, fault tolerant 
quantum computing, and high-impedance isolation for electrical current standards based on Bloch 
oscillations. 

PACS numbers: 74.50.+r, 74.81. Fa, 85.25.Am 
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Superinductors are desired for the implementation of 
the electrical current standards based on Bloch oscil- 
lations [1, 2], protection of Josephson qubits from the 
charge noise [3, 4], and fault tolerant quantum compu- 
tation [5, 6]. The realization of superinductors poses 
a challenge. Indeed, the "geometrical" inductance of a 
wire loop is accompanied by a sizable parasitic capaci- 
tance, and the loop impedance Z does not exceed aRq 
7], where a = 1/137 is the fine structure constant and 
Rq = h/4e 2 is the resistance quantum. This limita- 
tion does not apply to superconducting circuits whose 
kinetic inductance Lk is associated with the inertia of 
the Cooper pair condensate [8]. 

The kinetic inductance of a Josephson junction, Lk = 
($o/27r) 2 / Ej, scales inversely with its Josephson en- 
ergy Ej [8] ($ = h/2e is the flux quantum). The ki- 
netic inductance can be increased by reducing the in- 
plane junction dimensions and, thus, Ej. However, 
this resource is limited: with shrinking the junction 
size, the charging energy Eq — e 2 /2C (C is the junc- 
tion capacitance) increases and the phase-slip rate oc 
exp[— (SEj/Ec) 1 / 2 ] [8] grows exponentially, which leads 
to decoherence. Small Josephson energy (i.e., large ki- 
netic inductance) can be realized in chains of dc SQUIDs 
frustrated by the magnetic field [9, 10]. However, the 
phase-slip rate increases greatly with frustration, and the 
chains do not provide good isolation from the environ- 
ment. For the linear chains of Josephson junctions with 
Ej/Ec ^ 1, relatively large values of Lk (up to 300 nH 
['!]) have been realized in the phase-slip- free regime; fur- 
ther increase of the impedance of these chains is hindered 
by the growth of their parasitic capacitance. Also, the 
linear chains, as well as the nanoscale superconducting 
wires with a kinetic inductance of ~ 10nH//im [11, 12], 
are essentially linear elements whose inductance is not 
readily tunable (unless large currents are applied). 

We propose a novel superinductor design that has sev- 
eral interesting features. This circuit can be continuously 
tuned by a weak magnetic field between the regimes char- 



acterized by a low linear inductance and a very large 
nonlinear inductance. Importantly, the large impedance 
Z 3> Rq is realized when the decoherence processes as- 
sociated with phase slips are fully suppressed. This com- 
bination of strong nonlinearity and low decoherence rate 
is an asset for the development of high-performance su- 
perconducting qubits and controllable coupling between 
qubits. 

The studied circuit is a "ladder" of nanoscale 
Josephson junctions frustrated by the magnetic flux $ 
[Fig. 1(a)]. Each unit cell of the ladder represents 
an asymmetric dc-SQUID with a single "small" junc- 
tion with the Josephson energy Ejs in one arm and 
three "large" Josephson junctions with the Josephson 
energy Ljl in the other arm (the in-plane dimensions 
of both types of tunnel junctions do not exceed 0.3 x 
0.3 ^tm 2 ). The adjacent cells arc coupled via one large 
junction; the Josephson energy of the system, Ej(<p), re- 
mains an even function of the phase difference tp across 
the ladder at any value of the flux (the benefits of this 
symmetry are discussed below). 

For the Josephson junctions in the ladder "backbone," 
the Josephson energy is two orders of magnitude larger 
than the charging energy (the junction parameters are 
summarized in Table I). In this case quantum fluctua- 
tions of the phase across individual junctions are small, 
and the dependence Ej(ip) [see Fig. 1(b)] can be ob- 
tained from the classical computation that minimizes Ej 
with respect to the phases of superconducting islands 
at a fixed tp. At zero frustration ($ = 0) the energy 
Ej((p) is approximately parabolic over the relevant range 
of tp. With an increase of the magnetic field, the cur- 
vature of Ej((f) near tp = decreases. Provided the 
ratio of the Josephson energies for large and small junc- 
tions, r = Eji/Ejs, is not too large, the curvature 
d 2 Ej(tp)/dtp 2 \ v= o vanishes at some critical flux $ c . At 
this frustration the potential is strongly anharmonic (ap- 
proximately quartic) and the kinetic inductance, which 
is inversely proportional to the curvature of the potential 
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FIG. 1: (color online) Superinductor unit cells, potential energy, and phase diagram, (a) The unit cells of the tested device 
include "small" and "large" Josephson junctions with Josephson energies -Bjs and -Ejl, respectively. The "backbone" of the 
superinductor is shown as a bold line; the adjacent cells share large Josephson junctions. The cells are threaded by the same 
magnetic flux $; the phase difference across the device is <pi — (pi — (p. (b) The Josephson energy Ej(ip) of a six-cell ladder 
calculated within the quasiclassical approximation at $ = (dashed curve) and $ = $o/2 (solid and dash-dotted curves); 
Ejs = 3.5 K for all three curves, Eji, — 16.8 K for dashed and solid curves and -Ejl = 14 K for the dash-dotted curve. For the 
optimal value of the ratio Ejl/Ejs, denoted as r , the dependence Ej(<p) becomes flat near tp = at $ = $o/2 (solid curve). 
For smaller values of -Ejl/-Ejs, a double well potential is realized near full frustration, (c) "Phase diagram" of the ladders on 
the r-$ plane. The values of r are 4.1 and 4.5 for the ladders with TV = 6 and 24, respectively. The values of r — r for the 
studied devices are shown as horizontal lines. 



Lk oc [d 2 Ej(ip)/dip 2 ] 1 [8] in its minimum, diverges. 

The optimal regime of operation for the studied su- 
perinductor and supcrinductor-based qubits is realized 
when the curvature d 2 Ej(ip)/dip 2 \ tp= o vanishes exactly 
at full frustration ($ = $o/2). Indeed, because the en- 
ergy is an even function of the flux at full frustration, the 
device in this regime is insensitive in the first order to the 
flux noise [13]. This allows for simultaneous realization 
of the maximum inductance (i.e., the maximum fluctu- 
ations of the phase across the ladder) and the minimal 
coupling to the flux noise. This regime corresponds to 
the optimal value of the ratio Ejl/Ejs, which we denote 
as r Q . The r value, calculated within the classical ap- 
proximation (Eq = 0), varies between 3 for a single unit 
cell and 5 for a very long ladder. Quantum fluctuations 
renormalize Ej((p) and reduce r Q by <5r; for our devices 
Sr ~ 0.5-0.7. 

We have experimentally studied two types of ladders 
with the number of unit cells N = 6 and 24. Short lad- 
ders with N = 6 allow for direct comparison between 
experimental data and simulations based on the numer- 
ical diagonalization of the system Hamiltonian. Because 
these simulations are infcasiblc for a larger number of 
unit cells, the ladders with N = 24 were treated within 
the quasiclassical approximation. These longer ladders 
demonstrate the potential of our novel approach: their 
microwave impedance exceeds i?q by an order of magni- 
tude. 

The ladders and the readout circuits were fabricated 
using multiangle electron-beam deposition of Al films 
through a liftoff mask (see Supplemental Material IA for 
details). The in-plane dimensions of small and large junc- 
tions were 0.16 x 0.16 /im 2 and 0.30 x 0.30 /im 2 , respec- 



tively. The unit cell size was 3x5 /j,m 2 , and the flux 
$ = $0/2 was realized in the magnetic field B w 0.7 G. 
The ladders parameters are listed in Table I. Several de- 
vices with systematically varied values of r were fabri- 
cated on the same chip and inductively coupled to the 
same microstrip line. The devices could be individually 
addressed due to their different resonance frequencies. 

The effective kinetic inductance of the ladder, Lk = 
(u^Ck)^ 1 , was calculated from the measured fre- 
quency wrji of the |0) O |1) transition in the resonance 




FIG. 2: (color online) Schematic description of the on-chip 
circuit, (a) Topology of two unit cells of the ladder. Joseph- 
son contacts are formed at the intersections of the bottom 
(horizontal) and top (vertical) electrodes, (b) The circuit di- 
agram of the superinductor and LC resonators coupled via 
the kinetic inductance Lc of a narrow superconducting wire. 
The micrograph shows a ladder with six unit cells, (c) The 
on-chip circuit layout of two resonators inductively coupled 
to the microwave (MW) feedline. 
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TABLE I: Josephson junction parameters. Ecs and Ecl are the charging energies for the small and large junctions. 



Device 


Number 


Ejs, 


Ecs, 


Ejl, 




r = 


Ck, 


La, 


M* = o), 


= $o/2), 




of unit cells 


K 


K 


K 


K 


Ejh/Ejs 


fF 


nH 


nH 


nH 


1 


6 


3.2 


0.46 


14.5 


0.15 


4.5 


18 


0.4 


3.7 


130 


2 


24 


3.15 


0.46 


14.5 


0.15 


4.6 


5 


0.8 


16 


3 000 



circuit formed by the ladder and the interdigitated capac- 
itor. The capacitance Ck is larger than the capacitance 
of the interdigitated capacitor due to the parasitic capac- 
itance to the ground. The resonance frequency of this cir- 
cuit, referred below as the superinductor resonator, varies 
with the magnetic field by an order of magnitude (see 
below), whereas the bandwidth of the cryogenic pream- 
plifier and cold circulators in our measuring setup is lim- 
ited to the range of 3-10 GHz (the microwave setup is de- 
scribed in Supplemental Material IB). To overcome this 
limitation, the superinductor resonator was coupled to 
a linear lumped-element LC resonator with a resonance 
frequency ulc/^tt ~ 7 GHz via the coupling inductor 
Lq [a narrow superconducting wire, Figs. 2(b) and 2(c)]. 
The LC resonator is formed by an inductor (meandered 
2-/mi-wide Al wire) with L = 5 nH and an interdigitated 
capacitor (2-/im-wide fingers with 2 /im spacing between 
them) with C = 100 fF [Fig. 2(c)]. Both the superin- 
ductor and LC resonators are inductively coupled to a 2- 
port Al microstrip feedline with a 50 Q wave impedance. 
When the superinductor resonator is excited by a second- 
tone microwave frequency 0J2 , its impedance changes due 
to nonlinearity; this results in a shift of the resonance 
of the LC resonator probed by the first-tone microwaves 
with u>i w u>lc- The maximum shift occurs at the fre- 
quency 0J2 = ujqi that corresponds to the |0) |1) tran- 
sition. 

Figure 3 shows the resonance modes corresponding to 
the transition |0) <H- |1) in the superinductor resonators 
with the 6-cell and 24-cell ladders. The insets in Fig. 3 
show the avoided crossing between the lowest modes of 
the superinductor and LC resonators observed in the 
first-tone measurements (these avoided crossings illus- 
trate the strength of coupling between these resonators). 
The measured values of the lowest-mode frequency woi 
for the 6-ccll ladder arc in excellent agreement with simu- 
lations based on the numerical diagonalization of the cir- 
cuit Hamiltonian (see Supplemental Material IIB). The 
only fitting parameter in these calculations was the ra- 
tio r = Ejl/Ejs, which is within 15% of the designed 
value. The nominal junction parameters for the 24-cell 
ladder also agree with the quasiclassical simulation of the 
dependence cjoi(3>)- 

Superinductor applications require that the fre- 
quency w i, which corresponds to the "global" mode of 
the superinductor, should be much smaller than the fre- 
quency of its internal excitations, Wi n t. Far away from full 
frustration, the internal modes correspond to very high 



frequencies of an order of the Josephson plasma frequency 
(~ 100 GHz for our samples). At full frustration, Wj n t de- 
creases, but, according to our estimate, remains above 
10 GHz for the 24-cell device. 

The inductance at full frustration increases as r ap- 
proaching its optimal value r , which depends on the 
number of unit cells (see Supplemental Material IIB): 
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FIG. 3: (color online) Spectroscopic data for the ladders 
with N = 6 (a) and N = 24 (b). The resonance frequen- 
cies woi (<&)/2tt of the |0) o |1) transition measured in the 
second-tone experiments are shown by red dots. The hori- 
zontal dashed lines correspond to the resonance frequency of 
the LC resonators. The blue curves represent the fits based 
on numerical diagonalization of the circuit Hamiltonian for 
device 1 (a) and the quasiclassical modeling for device 2 (a); 
for the simulation parameters, see Table I. The gray scale in- 
sets show the microwave amplitude 1521 1 versus the first-tone 
microwave frequency uji/2iv and the normalized flux $/$o 
measured near the avoided level crossings. 
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r Q = 4.1 for N = 6 and r D = 4.5 for N = 24. Even 
though both devices 1 and 2 have similar values of r (i.e. 
nominally identical junction parameters), their proximity 
to the critical point is significantly different. For device 1 
with r — r Q m 0.4, the inductance at full frustration ex- 
ceeds that at zero field by a factor of 35. For device 2 
with r — r m 0.1, this increase exceeds two orders of 
magnitude, and the inductance at full frustration is 3 fill 
(for comparison, this is the inductance of a 3-meter-long 
wire). The total capacitance of the superinductor res- 
onator, Ck = 5fF (see Table I), includes the capacitance 
of the interdigitated capacitor (2 fF) and the parasitic 
capacitance of all wires and the superinductor to the 
ground (3fF) obtained by circuit modeling. By "mean- 
dering" the ladder and shrinking the tested moderately- 
sized (3x5 /im 2 ) unit cell, the parasitic capacitance can 
be reduced down to ~ 1 fF ; the impedance of such a 
ladder approaches 50kil at oj/2it = 3 GHz. 

The nonlinearity of the superinductor-based qubit in- 
creases dramatically with approaching the optimal work- 
ing regime (r = r D and <& = $o/2)- According to our 
quantum simulations, the nonlinearity factor 7 = (u>i2 ~ 
woi)/woi (uji2 is the frequency of the |1) <-t |2) transi- 
tion) for the 24-cell ladder (device 2) with r — r w 0.1 
approaches 40 % at full frustration. Strong, tunable quar- 
tic nonlinearity of the studied superinductor is an asset 
for the qubit design. In particular, strong anharmonic- 
ity enables fast qubit operations and qubit readout due 
to a large energy difference between the |0) -<->• |1) and 
|1) <-> 1 2) transitions [14]. 

In order to demonstrate the high quality of our su- 
perinductor, we have measured Rabi oscillations in the 
qubit formed by the superinductor and the capacitor Ck 
[Fig. 4(a)]. In these measurements, we have monitored 
the phase shift of the LC resonance while the superin- 
ductor resonator was excited by pulsed microwaves with 
LU2 = loqi- The quantum nature of these oscillations 
was verified by observing the linear dependence of the 
Rabi frequency on the amplitude of the microwave field 
[Fig. 4(b)]. The observed decay time of Rabi oscillations 
exceeded 1 /is and was limited by the energy relaxation 
time [cf. Fig. 4(c)]. The dominant source of energy re- 
laxation is the intentional inductive coupling to the LC 
resonator and the microwave fecdlinc. 

The intrinsic decoherence rate of the qubit is expected 
to be very low. Because the curvature d 2 Ej(ip)/d(p 2 
(which controls the position of energy levels) has a min- 
imum at full frustration, one expects that the flux noise 
does not affect the qubit in the linear order. Another 
common source of dephasing in a chain of superconduct- 
ing islands coupled by Josephson junctions is the phase- 
slip processes in combination with ubiquitous fluctua- 
tions of offset charges on the islands [15]. Because of 
the Aharonov-Casher effect [3, 16-18], these two factors 
result in the decoherence which is directly proportional to 
the phase-slip rate. In the studied devices this rate, being 
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FIG. 4: (color online) (a) Rabi oscillations of the population 
of the first excited level of the superinductor resonator in de- 
vice 1. The phase shift of the LC resonance was measured 
while the superinductor resonator was excited by the second- 
tone pulsed microwaves with U12 = woi- The data are shown 
for the phase ip = 27r<l?/$o = 0.947T. The solid line repre- 
sents the fit with the Rabi decay time 1.4 fis. (b) Dependence 
of the Rabi frequency on the amplitude of microwaves with 
u)2 = W01 . (c) The response of the LC resonator (measured at 
uji = ujlc) to the excitation of the superinductor resonator 
by a 0.4 /js second-tone pulse. 

proportional to exp[— ^Ejl/Ecl) 1 ^ 2 }, where c 2.5-2.8, 
is expected to be negligible due to a large ratio Ejl/Ecl 
(sa 100) for the junctions that form the ladder backbone 
[Fig. 1(a)]. Recent work [19] demonstrated that linear 
chains of Josephson junctions with Ej/Ec ~ 100 are 
phase-slip-free and exhibit inductances up to 0.3 /iH. 

We envision many applications for the designed su- 
perinductor: this element has the potential to reduce 
the charge noise sensitivity of Josephson qubits, enable 
implementation of the fault tolerant qubits, and pro- 
vide sufficient isolation for the electrical current stan- 
dards based on Bloch oscillations. The ability to trans- 
form this element from an inductor with an almost lin- 
ear response into a very nonlinear superinductor by tun- 
ing the magnetic field can facilitate controllable coupling 
between qubits. In a moderately or strongly nonlinear 
regimes, the superinductor-based resonator can also op- 
erate as a qubit with a low decoherence rate. Being com- 
bined with a small Josephson junction, the superinduc- 
tor can be used as an adiabatic switch — an element 
whose impedance changes exponentially with magnetic 
field [5, 20], which is crucial for the fault tolerant qubit 
operations [6]. 

We would like to thank V. Manucharyan and 
A. Zamolodchikov for helpful discussions. The work was 
supported by DARPA (HR0011-09-1-0009), NSF (DMR 
1006265), and ARO (W911NF-09-1-0395). 
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Appendix A: Experimental details 
1. Device fabrication 

The superinductor, the lumped-element LC resonator, 
and the microstrip fecdline line were fabricated within 
the same vacuum cycle using multi-angle electron-beam 
deposition of Al films through a nanoscale lift-off mask. 
To minimize the spread of the junction parameters, we 
have used the so-called "Manhattan-pattern" bi-laycr 
lift-off mask formed by a 400-nm-thick e-bcam resist (the 
top layer) and 50-nm-thick copolymer (the bottom layer) . 



In this technique, the circuit is formed by aluminum 
strips of a well-controlled width intersecting at right an- 
gles (thus, the "streets" and "avenues"). After depositing 
the photoresist on an undoped Si substrate and expos- 
ing the pattern with e-beam lithography, the sample was 
placed in an ozone asher to remove any traces of the 
photoresist residue. This step is crucial for reducing the 
spread in junction parameters. 

The substrate was then placed in an oil-free high- 
vacuum chamber with a base pressure of 5 x 10 -9 mbar. 
The rotatable substrate holder is positioned at an angle 
of 45° with respect to the direction of e-gun deposition of 
Al. During the first Al deposition, the substrate was posi- 
tioned such that the "streets" were orthogonal to the de- 
position direction. If the width of lines in the e-beam re- 
sist forming the "streets" is less than the mask thickness 
(0.45 imo), no metal is deposited in the "streets" during 
this deposition, whereas the "avenues" are fully covered 
with Al. Alternatively, if the lines forming the "streets" 
are wider than 0.5 /mi, both the "streets" and "avenues" 
are covered with metal in the first deposition. The lines 
narrower than 0.5 fim were used to form Josephson junc- 
tions, the wider lines were used for the fabrication of 
meander-shaped inductors. Without breaking vacuum, 
the surface of the bottom Al layer with a thickness of 
20 nm was oxidized at ~ lOOmTorr of dry oxygen for 5 
minutes. 

After evacuating oxygen, the substrate holder was ro- 
tated by 90° , and the second 60-nm-thick Al film was de- 
posited; this layer forms top electrodes of Josephson junc- 
tions, and increases the total thickness of the meander- 
shaped inductors and intcrdigitatcd capacitors. Finally, 
the sample was removed from the vacuum chamber and 
the lift-off mask was dissolved in the resist remover. The 
spread of the resistances for the nominally identical JJs 
with an area of 0.15 x 0.15 /im 2 did not exceed 10%. 



2. Measurement technique 

The microwave response of the coupled superinductor 
and LC resonators was probed by measuring both the 
phase and amplitude of the microwaves traveling along a 
microstrip feedline coupled to the resonators. A simpli- 
fied schematic of the microwave circuit is shown in Fig. 5. 
The cold attenuators and low-pass filters in the input mi- 
crowave line prevented leakage of thermal radiation into 
the resonator. On the output line, two cryogenic Pamtech 
isolators (~ 18 dB isolation between 3 and 12 GHz) an- 
chored to the mixing chamber were used to attenuate the 
~ 5K noise from the cryogenic amplifier. 

The probe microwaves at frequency u\, generated by 
a microwave synthesizer (Anritsu MG3694B), were cou- 
pled to the cryostat input line through a 16 dB coupler. 
These microwaves, after passing the sample, were ampli- 
fied by a cryogenic HEMT amplifier (Caltech CITCRYO 



6 




16 dB Coupler 



FIG. 5: Simplified circuit diagram of the measurement setup. 
The microwaves at the probe frequency Ui are transmitted 
through the microstrip line coupled to the superinductor and 
LC resonators. The microwaves at ui r , after mixing with the 
microwaves at ui, provided the reference phase 4>o- The signal 
at Wi is amplified, mixed down to an intermediate frequency 
uji —CJ r by mixer M2, and digitized by a fast digitizer (ADC). 
The second channel of the ADC is used to digitize the signal 
from mixer M2. The microwaves at the second-tone frequency 
u>2 were coupled to the transmission line via a 16 dB coupler. 



1-12, 35 dB gain between 1 and 12 GHz) and two 30 dB 
room-temperature amplifiers. The amplified signal was 
mixed by mixer Mi with the local oscillator signal at 
frequency w r , generated by another synthesizer (Giga- 
tronics 910). The intermediate-frequency signal A(t) = 
Asin(ilt + <po) + ^4n(*) (where A-^(t) is the noise term) 
at £1 = (uji —ijj r )/2it = 30 MHz was digitized by a 1 GS/s 
digitizing card (AlazarTech ATS9870). The signal was 
digitally multiplied by sin(f2t) and cos(r2f) averaged over 
~ 10 6 periods, and its amplitude A (proportional to 
the microwave amplitude IS21I) and phase <fi were ex- 
tracted as A = ([A(t) sin(fit)] 2 + [A(t) cos(f2i)] 2 ) 1/2 and 
<t> = &i'ctaiL{([A{t)sin{nt)} 2 )/([A{t)cos{flt)} 2 )}, respec- 
tively. The reference phase </>o (which randomly changes 
when both kq and u>2 are varied in measurements) was 
found using similar processing of the low-noise signal pro- 
vided by mixer M2 and digitized by the second channel 
of the ADC. This setup enables accurate measurements 
of small changes 4>~4>o unaffected by the phase jitter be- 



tween the two synthesizers. The low noise of this setup 
allowed us to perform measurements at a microwave 
excitation level of — 133 dBm which corresponded to a 
sub-singlc-photon population of the tank circuit. In the 
second-tone measurements, the superinductor resonator 
was excited by the microwaves at frequency 0J2 propagat- 
ing along the same microwave feedline that was used for 
the microwave transmission at uq. 



Appendix B: Theoretical analysis 

Here wc explain in detail the nature of the large in- 
ductances of the Joscphson ladders studied in this work, 
the role played by the field induced frustration and the 
numerical computations which results we compared with 
the data in the main text. We begin our discussion with 
the analysis of the classical (the Coulomb energy Ec = 0) 
model. We then discuss the computation that takes into 
account quantum fluctuations induced by a non-zero Eq ■ 



1. Classical analysis 

The main idea of the superinductor design becomes 
more transparent in the case of a fully frustrated ($ = 
<&o/2) long ladder shown in Fig. 6(a). As wc explain 
below, the fully frustrated regime is also very important 
experimentally because in this regime the quantum states 
of the ladder are less sensitive to the flux noise. 

It is convenient to choose the gauge in which the mag- 
netic field induces the phases 27r<I>/$o on smaller junc- 
tions. Because at full frustration these phases are equal 
to 7r, the shift of the ladder by one unit cell combined 
with its rotation about its axis transforms the ladder into 
itself, so the problem becomes translationally invariant. 
Furthermore, a large length of the ladder allows to ne- 
glect the boundary effects and focus on translationally 
invariant solutions. 

Consider a small portion of the fully frustrated ladder 
shown in Fig. 6(a). The translational invariance im- 
plies that the solution is described by two phases across 
large junctions: a at the "vertical" junctions and j3 at 
the "horizontal" ones. The total phase across the whole 
ladder is tp — N(a + /3), where N is the number of rungs. 
The energy per unit cell that contains two large and one 
small junctions is given by 

E(a,/3) = --Ej L (cosa + cos/3) — Ej S cos(7r — 2a — /?). 

(Bl) 

The dependence E(a, j3) is shown in Fig. 7(a). At Ejj, » 
Ejs the first term in Eq. (Bl) dominates and E(a, (3) has 
a minimum at (a,/3) = (0,0). The expansion near this 
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FIG. 6: (a) A small portion of a long fully frustrated Joseph- 
son ladder. Large junctions with the Josephson energy Ejl 
and the Coulomb energy Ecl -C -Bjl are shown in red 
thick lines; they form the ladder backbone. Small junc- 
tions with the Josephson energy _Bjs and the Coulomb energy 
Ecs "C Ejs are shown in cyan thin lines. The frustration of 
each cell is controlled by the dimensionless parameter &/<&q. 
(b) Schematic diagram of a six cell ladder (device 1). The 
ends of the ladder are grounded: the phase difference across 
the ladder is ip. 



point gives 



Ejl — 4-Ejs — 2-Ejs 




a 


—2Ejs -E.il — Ejs 




(3 



(B2) 

which is a quadratic form with eigenvalues Ejl and 
-Ejl — 5i?jg. The second eigenvalue changes sign as the 
ratio r = Ej L /Ej S decreases, signaling the instability 
and appearance of a nontrivial ground state. The crit- 
ical value of the ratio, r = 5, corresponds to the op- 
timal working point, as explained below. At this point 
the function E(a,f3) becomes flat along the eigenvector 
(2,1). We introduce new coordinates in the "flat" and 
"steep" directions: 7 = (2a + 0)/VE, 5 = (a - 20) /VE. 
Figure 7(b) shows the energy plotted along the flat di- 
rection. Neglecting the phase deviations in the steep di- 
rection, the total phase is related to 7 by <p = 3N/^/Ej. 

At r > r Q the ground state values of phases 7 and ip, 
70, and (po, are zero. Deviations of the total phase from 
tpa = result in the quadratic increase of the energy 



18N 



At the optimal point (r = r a ) these deviations van- 
ish, and the inductance, which is inversely proportinal 
to d 2 E/d(p 2 , becomes infinite in the quadratic approxi- 
amtion and is defined by the non- vanishing quartic term 



E^(cp) 



625 - 17r 
48600AT 3 



EjsV 



90N 3 



EjsV ■ 



At smaller r < r Q the ground state corresponds to a non- 
zero value of |7o(^)| oc \/r a — r that translates into a non- 
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FIG. 7: (a) Energy per elementary cell E(a, 0) (Bl) as a func- 
tion of a and /3 at the optimal point r — r a . The black line 
corresponds to the flat direction parametrized by 7. (b) The 
dependence of energy in the flat direction for different values 
of r. (c) Phase diagram of the model in classical approxi- 
mation in dimensionless variables r = Ejl/Ejs and &/<&o- 
Quantum fluctuations discussed in section B 2 shift the tran- 
sition line down but affect neither the long range order nor 
the phase transition. 



zero <fo = ±3Nj /VE. The appearance of a non-zero 70 
is equivalent to a phase transition. 

The translational invariance discussed above is vio- 
lated away from full frustration. At arbitrary frustra- 
tion, the ladder remains symmetric under the translation 
by two unit cells. Denoting the phases across the large 
junctions in the first cell by a, and in the second cell 
by a',/3' we obtain the classical energy per two rungs 

= -Ejl [cos a + cos (3 + cos a' + cos f3'\ 



-Ejs 



cos( 



2tt$ 
3>n 



2tt$ 



cv+ - 13) + cos(— — +a + + f3') 



( I> ( ] 



(B3) 



where a + = a + a' . This energy does not change un- 
der the combined transformation a — > a', /3 0', and 
$ — > — $, which correspond to the shift by one unit 
cell and simultaneous change in the sign of the mag- 
netic field. The total phase difference along the ladder 
ip = N(a + +[3 + )/2 is invariant under this transformation. 
Time inversion invariance implies that the ground state 
energy does not change under the simultaneous change in 
sign in both $ and tp, the symmetry under the combined 
transformation implies that the ground state energy is an 
even function of ip at any field. This conclusion is based 
only on the symmetry of the problem and thus remains 
valid when quantum fluctuations arc taken into accound 
(see section B 2). 
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The energy (Bl) is minimal at a_ = a' — a — and 
/3_ = (/?' — /3)/2 7^ that has to be determined from the 
minimization of (B3). The minimum in cv_./3_ is always 
steep, so fluctuations in this direction can be neglected. 
In this approximation the general expression for the en- 
ergy (B3) is reduced to the form (Bl) with the effective 

#js ff) =£j S cos(27nr + /?_), 



where x = $/$o - 1/2 



At r > r Q the ratio Ej-l/Ej ^ that controls the insta- 
bility is large: Ej^/Ejg > r > r a , so the ground state 
corresponds to zero phases. At r < r a the instability oc- 
curs at non-zero x c , defined by Ej-^/Ej^ (x c ) — r Q [see 
Fig. 7(c)]. At the point of instability we can determine 
the value of assuming that a+ = /3+ = 0, we get 



P—(x) — arctan 



sin(27ra;) 



r — cos(27rx) 



(B4) 



The effective reduction of Ej$ implies that at r > r Q 
the transition from the single minimum at 7 = to dou- 
ble minimum at ±70 occurs as a function of the field at x c 
where 



cos[27rx c + /3_(a; c )] 



(B5) 



The transition line defined by Eqs. (B4) and (B5) sepa- 
rates a symmetry-broken (ordered) phase from the phase 
with unbroken symmetry. We show the phase diagram 
obtained by the numerical solution of these equations in 
Fig. 7(c). The equation for the instability can be solved 
analytically for r Q — r <C r a (i.e. x <C 1): 

r ( r 

x c = r arccos — 

27r(r + l) \r Q 

At the critical flux $ c = (1/2 ± x c ) $0 the quadratic part 
of the energy vanishes which implies infinite inductance. 
Thus, the system can be tuned to the infinite inductance 
either by realizing r = r Q at full frustration, or by ap- 
proaching <f> c at r < r a . 

In the absence of noises and disorder, these two meth- 
ods of realization of the infinite inductance (and the 
phase transition) are equivalent. However, the real sys- 
tems are subject to flux noise. Because the energy of the 
ground state is quadratic in flux deviations (5$, the effect 
of the flux noise on the system is minimized if the infi- 
nite inductance is realized at full frustration, or close to 
it. Thus, the optimal working point corresponds to the 
ratio Ej L /E JS = r Q . 

The finite length of the ladder brings in two effects 
associated with the boundaries. First, this may violate 
the symmetry which resulted in E(cp) = E(—tp) for all 
frustrations; this is crucial for the realization of infinite 
inductance at $ c ^ $0/2. However, this effect is ab- 
sent for the odd number of rungs (this was the reason 
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FIG. 8: (a) Inductance as a function of frustration (external 
magnetic field) in the classical approximation for a three-rung 
ladder at r — r a . (b) The inductance as a function of the ratio 
r = Ejl/Ejs for ladder with 2, 3, 5, and 7 rungs. Optimal 
ratios r : 3 for 2 rungs, 4 for 3 rungs, 4.61 for 5 rungs, and 
4.80 for 7 rungs. 



for selecting the even number of unit cells in the studied 
ladders). In this case the ladder remains symmetric un- 
der the rotation by it about its center. The rotation does 
not change the magnetic fluxes but changes ip — > —(p, 
so in this case the function E(ip) remains even for all 
fluxes and r. Second, the finite-size effects decrease the 
value of the critical r at which the instability occurs at 
full frustration. To compute them for finite ladders we 
numerically computed the energy E(ip) for the ladders 
containing from 2 to 7 rungs. The dependences of the 



for different 



ip=ipo 



inductance L K (<!>,r) = (d 2 E j dip 2 ) 

ladder lengths are shown in Fig. 8. The value of ipois de- 
fined by the minimum of the ground state energy E(tp) as 
a function of the phase ip. The inductance increases with 
frustration and diverges at full frustration $ = $o/2 for 
r = r Q . For the ladders with six unit cells, the finite size 
effects lead to a small decrease of the optimal value to 



2. Quantum analysis 

Non-zero (but small) charging energies result in quan- 
tum fluctuations of the phase of individual islands, Qi. 
These fluctuations have three effects. First, they might 
result in a phase slip in which the phase difference rjij = 
6i — 8j between two islands connected by a large Joseph- 
son junction changes by 2-7r due to a quantum tunneling 
process. Second, they may result in renormalization of 
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FIG. 9: Solid lines in panels (a) and (b) show the inductance 
computed in the classical approximation for a three-rung lad- 
der as a function of frustration at r = r (a) and as a function 
of r for $ = "3?o/2 (b). Dashed lines show the results of the full 
quantum computation for the same ladder. Quantum fluctu- 
ations shift position of the optimal point from its classical 
value (r = 4) to r » 3.8. Dashed line in panel (a) shows the 
result of the full quantum computation at r = r which takes 
into account the shift of r due to quantum fluctuations. 



the effective energy in the flat direction discussed in sec- 
tion B 1 and the shift of the phase boundary in Fig. 7(c). 
Finally, they may have a very large effect on the proper- 
ties of long ladders close to the transition point at which 
critical fluctuations are expected to become important. 

In the absence of frustration the amplitude of a 
phase slip across a single large junction, i? p h ~ 

3.8E c ^Ej(^ exp(— ■\/8£ , jl/£ , cl), is exponentially small 
for the experimental values of Ejl/Ecl ~ 100. The frus- 
tration decreases the tunnel barrier resulting in larger 
phase slip amplitudes. However, this effect is numeri- 
cally small. Namely, the dominant transitions correspond 
to the process in which the phase ft changes between 
and 2ir [sec Fig. 7(a)], the action for this process is 
S « 2.5^/EjiJEcl, Thus, the phase slip amplitudes 
remain negligible for the ladders with experimental pa- 
rameters. 

Small quantum fluctuations of the phase differences 
smear the classical potential. As a result the large induc- 
tance and the phase transition are realized at a smaller 
value of r . To estimate this effect we have diagonal- 
ized numerically the Hamiltonian describing small lad- 
ders. This diagonalization is convenient to perform in 
the charge basis in which the general Hamiltonian of the 
Josephson circuit acquires the form 



lA b/b • I I.e. ) 



(2e) 2 J2(C- 1 ) ij (n i + q i )(n j +q j ), (B6) 



where matrix Jy describes Josephson couplings between 
islands and Cij is the capacitance matrix of the sys- 
tem. For large junctions J»,j+i = J5jl, CVi+i = Cl, 
and Ecl = e 2 /2C*L; for small junctions JiA+3 = Ejs, 
Cm+3 = Cs, and Eqs = e 2 /2Cs. The operators h^ 1 
increase (decrease) the number of Cooper pairs on each 
island: = |n,- ± l)j(n,|. The junction matrix Jij is 
symmetric and its diagonal terms are zeros. The capac- 
itance matrix CV,- is positively defined with positive di- 
agonal elements Cn and negative off-diagonal elements 
Cij, i 7^ J, so that its inverse (C -1 )^ is positively de- 
fined and its elements are all positive. The phases Aij 
are induced by the magnetic field, the sum of ^4's over a 
closed loop I of n Josephson junctions corresponds to the 
flux penetrating this loop. 

For the numerical diagonalization we had to limit the 
number of charging states, n q at each island. For the pa- 
rameters similar to the experimental ones it is sufficient 
to keep n q = 12 to get the results with better than 1% 
accuracy. The effect of the quantum fluctuations on the 
inductance of a three rung chain with Ej-^/Eql = 40 is 
shown in Fig. 9. We observe that quantum fluctuations 
reduce r a by Sr m 0.2 down from its classical value, but 
the behavior of the inductance as a function of the frus- 
tration remains essentially the same. The importance of 
quantum fluctuations increases with the ladder length, so 
this shift becomes Sr w 0.5 for the ladders with six unit 
cells. 

The inductance L can be calculated as the inverse cur- 
vature of the ground state energy E(tp). Alternatively, 
it can be defined as L = l/ui^C for a superinductor 
shunted by a capacitor C. Here the frequency loqi is 
the excitation frequency of the superinductor resonator 
luqi = Ei — Eq. In the limit of a large capacitance C 
both methods give the same result (for a linear circuit, 
the results of both methods are identical for any C) . 

For the detailed comparison with the data for device 1 
we have diagonalizcd the Hamiltonian of the ladder con- 
nected to a large capacitor. In this Hamiltonian we 
used the values of the charging energy of small junctions 
and their Josephson energy extracted from the indepen- 
dent measurements of the Cooper pair transistors con- 
taining nominally identical junctions with Ecs = 0.46 K 
and Ejs = 3.2 K. The value of the charging energy of 
the large junction that was calculated from the ratio of 
the areas of small and large junctions: Ecl = 0.15K. 
We used the designed value of the large capacitance 
Ck = 18 fF. We then diagonalizcd the Hamiltonian 
of the ladders containing up to six rungs, determined 
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FIG. 10: Classical and quantum B(r). From right to left: the 
results of classical computation and from the full diagonaliza- 
tion of Hamiltonians for three and four rungs ladders. 

the excitation frequency of the supcrinductor resonator, 
woi($) = E\ — Eq, and extrapolated the result to 7 rungs 
that cannot be computed directly due to an enormous 
size of the Hilbcrt space. For the experimental values 
of the charging energies r Q is shifted down to approxi- 
mately 4.1. The 7- rung devices studied in this work (e.g. 
device 1) are characterized by the ratios r that are sig- 
nificantly above r OJ so the extrapolation to seven rungs 
works reasonably well. We found that in this parameter 
range the form of the dependence aj i( c E>) in the interval 
0.4 < $/$o < 0.5 is very similar for different r, whilst 
the absolute values of woi(0.4$ ) an d woi(0.5$) change 
dramatically as a function of r. This dependence is plot- 
ted in Fig. 3(a). Fitting the ratio w O i(0.4$ o )/woi(0.5$o) 
to experimental values for device 1 gives for the ratio 
r = 4.5 that is close to the designed value of 4.0. 



The energy of the first excited states of longer device 2 
was obtained in a two step procedure. In the first step we 
used exact diagonalization to find the energy functional 
of small quantum ladders (with 3-5 rungs) forming the 
loop penertrated by flux $ = 4>o/2 and fit it to the gen- 
eral form expected in vicinity of the optimal point: 

E{<p) = W + %r W- (B7) 

Here B q (r) and C q (r) are numerical coefficients that dif- 
fer somewhat from their classical values, B c \ and C c \. 
Comparison of the results of the computation with clas- 
sical result shows that the main effect of quantum com- 
putation is a shift of the critical value of r (see Fig. 10), 
Bq(r) « B c \(r + 8r) with 6r ps 0.45 and a modest renor- 
malization of C q (r) sa 0.7C c i(r = r Q ) that can be approx- 
imated by the constant in the relevant regime. 

The procedure outlined above allowed us to deter- 
mine the renormalization of the classical energy by short 
scale quantum fluctuations. In the second step we used 
Eq. (B7) to compute the potential energy of the long lad- 
der of device 2. Solving finally the quantum problem of 
the phase fluctuations across the whole ladder character- 
ized by the Hamiltonian 

H = E{lp) + 4E c q 2 
we find the low energy excitations of the device. The 
results are shown in Fig. 3(b). 



